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The projector onto single quantum map eigenstates is written only in terms of powers of the 
evolution operator, up to half the Heisenberg time, and its traces. These powers are semiclassically 
approximated, by a complex generating function, giving the Husimi distribution of the eigenstates. 
The results are tested on the Cat and Baker maps. 

I. INTRODUCTION 

The relation between the classical invariant structures of a dynamical system and its quantum counterparts, the 
eigenenergies and eigenstates, is not completely understood and is an intense subject of investigation. When the 
system is integrable the eigenenergies are determined by the well-known Bohr-Sommerfeld or EBK quantization rules 
and the eigenfunctions have a very intuitive phase space structure: they are essentially localized on the quantized 
tori and they present characteristic quantum features -interference fringes for the Wigner function or zeros for 
the Husimi function away from them. No such clear picture is available for the eigenstates of chaotic systems 
and much work |^-^ has gone into trying to disentangle the purely statistical aspects of the eigenfunctions from the 
specific dynamical features such as scars and localization. 

Although a lot of recent work deals with the study of the statistical aspects of chaotic wavefunctions 0| we concen- 
trate here on the purely dynamical aspects and explore how well the semiclassical formulae can reproduce the details 
of individual eigenfunctions. 

With regards to spectral properties the theoretical tool for this analysis is the Gutzwiller trace formula |^ relating 
the spectral density to the periodic orbits. As is well known, this relationship is fraught with uncertainties both due to 
the approximations involved, as to the fact that it diverges where it is supposed to locate the single eigenenergies and 
to the sheer complexity of adding the contributions of an exponentially increasing number of periodic orbits. Some 
of these difhculties can be overcome by resummation techniques |^ that limit the exponential increase by partially 
restoring the unitarity of the propagator, lost in the semiclassical approximation. 

In the present work we extend these resummation methods to the calculation of eigenfunctions of chaotic maps. 
We utilize Fredholm theory to construct a projector on single eigenstates whose expression is directly amenable 
to semiclassical approximations. These approximations share the same advantages - and the drawbacks- of the 
resummation techniques available for spectral properties. 

In section 2 we construct explicitly the projector on single eigenstates. This is achieved exactly in terms of a finite 
sum on powers of the quantum map propagator and its traces. The cut-off time is half the Heisenberg time due to the 
explicit imposition of unitarity. In section 3 these formulae are implemented semiclassically by a) approximating the 
traces in terms of actions and instabilities of periodic orbits in the usual way, and b) by computing the propagator in 
the coherent state representation in terms of specially devised complex generating functions. 

The final result is a semiclassical expression for the Husimi function of an eigenstate in terms of properties of a 
finite set of periodic orbits. 

In section 4 we apply the method to the cat map -where the semiclassical approximation is exact - and to the 
baker map where we can test the approximation. 

The formulae are here applied as a test to simple maps, but they can be easily adapted to more realistic Hamiltonian 
systems by selecting an appropriate Poincare section . 

II. PROJECTOR ON SINGLE EIGENSTATES 

For a map, the density of Floquet states can be expressed in the usual way in terms of traces by the formula 

d{e)^- + Y.[trU-e-^"^]+C.C. (1) 

n— 1 
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1/N plays the role of Plank's constant /i, n is the discrete time, and U is the unitary quantization of the map. The 
density of states results in a 27r periodic function with unit (5-spikes at the eigenangles efe, defined by 

f/l^fe) = e^'lV'fe) • (2) 

Semiclassical approximations using this formula require the calculation of traces of powers of the map for large 
values of n which are in turn related to classical orbits of the map by the Gutzwiller- Tabor Q trace formula 

bn = trir" ~^A^ exp [i{S^/h - /i^7r/2)] , (3) 

7 

where 7 label the n-periodic orbits of the map, 5*^ is the action, fi-y is the Maslov index, and is a coefficient related 
to the monodromy matrix. The approximation is valid for fixed n in the limit as iV — > 00 (but not vice- versa). 
A more efficient scheme starts from the spectral determinant 

P(s) = det(/ - sU) (4) 

which vanishes at s = Sfc = e~*'^'° and which can be expanded as 

JV 



P(,s) = ^/3„s" . (5) 

The coefficients /3„ can be computed recursively from the traces as 

1 " 

and therefore when computed semiclassically as in (3) they are expressed in terms of linear combinations of periodic 
orbits (pseudo orbits jT2|) of period up to n. 

The advantage over the calculation in term of traces is that only N coefficients need to be calculated. Moreover as 
a consequence of unitarity the symmetry 

f3N-j ^ {-i f detU 0j (7) 

can be imposed, so that only N/2 coefficients are needed and therefore only periodic trajectories for times up to N/2 
(i.e. half the Heisenberg time) are involved for a semiclassical calculation. This is now an optimal encoding of the 
eigenvalues in terms of traces, with the N real phases ek expressed exactly in terms of the N/2 complex coefficients 

fin 

It is convenient to define another function [ p^ , proportional to the spectral determinant, which has the same roots. 
Utilizing the symmetry (|^), for an even N, this function can be written 

When s is on the unit circle 1/s = s, and then the second sum is the complex conjugate of the first. Then the 
function Z{s) becomes real on |s| — 1. 

Unitarity also implies that the spectral determinant has all its zeros on the unit circle. However when the coefficients 
are computed semiclassically this is no longer guaranteed. If the symmetry (^ is imposed the resulting spectral 
determinant becomes self-inversive, which is a necessary-but by no means sufficient- condition for unitarity [ p^ . 
Self-inversive polynomials have the following properties [|l5| : 

1) All zeros are either on the unit circle or symmetric with respect to it. Thus as a function of a continuous 
parameter the zeros of a self-inversive polynomial can only leave the unit circle in pairs by degenerating on it. 

2) If P(s) has only single roots then dP/ds has no zeros on the unit circle. 
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A. Green Operator 



Very similar techniques can be implemented for the calculation of matrix elements and therefore of eigenfunctions. 
To obtain the eigenfunctions of the unitary matrix tj, we define a Green operator 

<}W = m. (9) 

I — sU 

which, as an analytic function of s, has poles at s = Sfc = e"*"^*, and the residues arc proportional to the projectors 
Pk = |V'fc)(V'fc|- To obtain normalized projectors, we instead define a normalized Green operator 

tr[G{s)] 

At the poles, the singularities cancel and we obtain g{sk) = Pk- ,f{U) is left arbitrary for the moment, and specific 
choices -made below- will be used to obtain formulae with difi'erent properties. 

We can expand the preceding expression utilizing the transpose cofactor matrix to write the Green operator 

9is) = ^(^j^'(^r-^j . (11) 

tr[fiU)Ct{I-sU)] 
The transpose cofactor matrix can be written as a finite power series in s 



C*(/-sC/)= ^,s"l„ , (12) 

n 

^„ = ^Af/"-V (13) 

i=0 

By the unitarity of U, the X„ matrices satisfy the following symmetry 

Xn-o = -U\-l)^ deiU X]_^ . (14) 

Utilizing this symmetry two choices of f{tf) lead to expansions with useful properties: 
a) For f{U) = U, for even A'', the corresponding Green function is 

. ^ t.:lv det u s^-^ y.:lv n .^g. 

E„To-^ tr (UX^) - det U E^iV tr (xl) ' 
Replacing the expansion of the X„ operators in powers of U, we obtain 



EfJt' Ciis) tr (U^+^) - det U s^-i Et^'o Ci(l/s) tr 



W 



in which the projector is expressed as a finite combination of the first N/2 powers of the propagator. 
The coefficients Cj are obtained from the /?„ 

N/2-1 

Ci{s) = J2 Pn-i . (17) 

n—i 

This normalized Green operator has no singularities as s goes around on the unit circle. This can be demonstrated 
because the denominator is equal to —dP/ds (see Appendix A). Then by the second property of the self-inversive 

polynomials it is guaranteed that this denominator never takes the zero value for s over the unit circle. This property, 
guaranteed by unitarity, will remain true -even when the traces and propagator are calculated semiclassically- as long 
as self inversiveness is maintained. 
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b) For f{U) = i + sU , G{s) in (9) is the Cayley transform of U and is therefore hermitian on the unit circle. In this 
case 



9b{s) 



s" Y 



with 



Yn — Xn + UXn-1 — 2X„ — /3„/ 



(18) 



(19) 



These new operators satisfy a simpler symmetry relation than the preceding X„, analogous to the symmetry of the 
Pn coefficients 



Yn-j = -{-if dct UyJ . 
In this case, for even N, the normalized Green operator can be written as 



9b{s) 





)-St(l/5) 


tr 




~tr 


St(l/s-)' 



(20) 



(21) 



where 



tis) 



(iN/2 



N/2-1 



\Yn,2 + E ^"""^'^ 



(22) 



In this expression, more symmetrical than (|l6|), the numerator is anti- hermitian on the unit circle, and the denomi- 
nator is purely imaginary. The basic operator S(s) is a Fourier transform of the Yn operators and their traces. When 
utilizing this finite representation, they then become the basic objects of study in the time domain, in place of the 
propagator and its traces. 

When f7 is a finite matrix, (jais) and gtis) are both rational analytic functions of s. They are identical at the 
eigenvalues s = Sk (on the unit circle). However their zero and singularity structure is quite different. 

In particular ga{s), as noted, has a denominator which does not vanish on the unit circle, even when the coefficients 
are computed approximately, as long as self-inversiveness is imposed on them. The operator gai-s) then, has no 
singularities in a small strip enclosing the unit circle and can be followed continuously from eigenvalue to eigenvalue. 
A disadvantage is that ga{s) is not hermitian (except at the values s — Sk)- On the contrary gb{s) is hermitian on 
the unit circle, but its singularities (the zeroes of the denominator) are precisely on the miit circle, in between two 
successive eigenvalues. 

It should be noticed that, although for simplicity all the above formulae are written for even N, they are easily 
adapted to the odd case, with minimal changes. 

The preceding operator expressions for gai-s) and gtis) provide a common structure that underlies all further 
semiclassical approximations. Depending on the representation used to evaluate them, they can yield results for 
transition probabilities, phase space distributions, response functions or wave function correlations. This is achieved 
by linking the spectral properties of the map to a finite number of traces and matrix elements of the propagator. 
To apply them one needs the semiclassical evaluation of the propagator traces as in (^)- which are representation 
independent- and also of the semiclassical propagator in the chosen representation. For example if they are computed 
in the coordinate representation they provide expressions for If the propagator is computed in the Weyl 

representation |l^ they yield the Wigner distribution of eigenf unctions jlj] . In all cases the unitarity of the propagator 
- which is generally only approximate in semiclassical treatments - is partially restored by the imposition of the 
symmetries (0,Q|o|). In the next section the above formalism is developed in the coherent state representation to 
obtain expressions for semiclassical Husimi distributions. 



III. SEMICLASSICAL APPROXIMATION IN THE COHERENT STATE REPRESENTATION 

The formulae derived in the previous section are exact and for maps can, at first sight, seem just a complicated 
way of looking at a simple algebraic eigenvalue problem. However, they are prepared in such a way as to make the 
semiclassical transition as clear as possible and to unveil the necessary approximations involved in such a transition. 
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The coherent state mean value of g{s) is 



7i(s, z, z) 



{z\9{s)\z) 
{z\z) 



(23) 



This is a phase space distribution, analytic in s, which becomes the Husimi distribution of the eigenstates \ipk) at the 
values Sfe = e""^*". 

Its properties, reflecting the zeroes and singularities of g{s) can be followed continuously in the s complex plane. 
Its calculation involves the N/2 coherent state return amplitudes 



C„(z,z) 



in addition of the N/2 traces 



{z\U"\z) 



bn — trU^^ — dz dz Cn{z, z) . 



(24) 



(25) 



The first admits immediately a well known semiclassical approximation in the time domain, namely the Gutzwiller- 
Tabor trace formula (|^) . The stationary phase approximations involved in the derivation are of course valid for fixed 
time n and for ?i ^ 0, or equivalently, N oo, and their use in the regime n « N/2 needed for the exact calculation 
is questionable. 

The semiclassical calculation of the return amplitudes Cn{z, z) can be done by stationary phase calculation of the 
integral 



where 



(2/k/>(9/|f^"l*>(ft|2j) % dqf , 



'ttTi 



is the coherent state in coordinate representation, and the propagator has the usual Van Vlcck form 



^J2-Kin ^ 



dqtdqf 



exp 



(26) 



(27) 



(28) 



The sum must be made over all the paths in the phase space of the corresponding classical map that connect qi with 
qf in n steps. Additional phases connecting the different branches, the Maslov indices, must be added to Sj"'\ 

The calculation of the integral ( p6| ) with the Van Vleck approximation ( p8| ) can be done as usual by stationary phase 
and the result has the structure of another Van Vleck formula with complex coordinates for a generating function. 



{zf\U-\z,)^Y.\ -'S- exp 



--Gz[z^,zf) 



dzidzf 

This new classical generating function is a complex Legendre transform of the action S 

1 i 
G3^~S+ -^{Pfqf-Piqi) + -^{q} + p} + q'^ + pf) , 

or equivalently, in the complex coordinates 

z = -L(,-.p) 
z = ^ (g + «p) , 



(29) 

(30) 

(31) 
(32) 



it can be written as 
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-S 



Zf + 2zfZf 



f> 4 

This function depends explicitly on Zi,Zf variables, and yields the equations of motion 



.9G3 

' dz. 



: dG3 

dzf 



(33) 



(34) 



These are algebraic equations that define the map implicitly provided 



^ where the sum must be made 



dzidzf 

over complex paths that connect the point Zi to Zf. These formulae show that the return amplitude C„ has as 
exponent the function 6*3(2, z) — izz which using (|3^ ) can be shown to be stationary at the periodic points. Quadratic 
approximation of the generating function at these points yields the coherent state return amplitudes as a sum of 
Gaussian peaks centered on them 



Cn{z, z) = exp 



lg^{zrz) + y. 



(35) 



1 [V^ 



(z — Z^Y + [u^ — 1)(^ ~ ^7)(^ " 



Z^) 



Z,f 



(36) 



where z^,z^ are the position of the periodic points, 5-y are the corresponding actions, and u-^^v^ and are the 
elements of the complex monodromy matrix defined by 

/ dzf/dzi dzf/dzi \ _ f u v 
\^ dzf/dzi dzf/dzi J ~ u 

calculated on the corresponding periodic point. 

Using (|3^) for the return amplitudes and (|^) for the traces, the Husimi distributions of eigenstates are explicitly 
written in terms of properties of periodic points of all periods up to half the Heisenberg time N/2. 

Although the semiclassical approximations are available in the time domain for the traces b„ and the return 
amplitudes C„(z, z) it is clear from our formulae that the essential ingredients for the calculations are the coefficients 
/3n and the operators X„ (or y„). These are related by (||) and (|l3|) to &„ and C„ and therefore can be recursively 
computed from them. This procedure yields the /3„ as linear combinations of products -the so called pseudo-orbits- of 
terms involving many periodic orbits of different periods and with signs implying many possible cancellations between 
long and short orbits. The same situation prevails for the Xn (or The expressions yield linear combinations of 
return amplitudes at many different times. In fact X„(z, z) give the phase space distribution of pseudo-orbits and it is 
easy to show that = {N — n) trXn ■ In view of the cancellations involved it would be very desirable to have a direct 
semiclassical derivation -both for /3„ and X„- that would not go through the indirect procedure of approximating 6„ 
and Un, and then recursively computing (3„ and X„. However to our knowledge such a derivation is not available. 




IV. NUMERICAL CALCULATIONS 



We will test the methods proposed on two simple well known maps: the cat map, and the baker map. 
For the cat map 



Pf ) \c d I \p. 



{mod 1) (38) 



with a, b, c and d integer numbers, and ad — be — 1, the Van Vleck formula in the qi, qj representation is exact p9| ] 
and the corresponding one in the coherent state representation can be computed explicitly. 
The transformation to complex coordinates leads to the map 



Zf 



U V \ I Zi 
V U I \ Zi 



(39) 



with 
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u=^[{a + d)+i{b-c)] , (40) 
v^l[ia-d)+zib + c)] . (41) 

and its quantization with ( p8|) is obtained exactly by Gaussian integration. The coherent state return amphtudes are 
computed with the formulae (|3^) and ( |3^ ) , with the advantage that the coefficients of the monodromy matrix are the 
same for all the periodic points of a given period. Since the cat map is defined on the torus we use periodized coherent 
states 1^. With this result we can test our formulae keeping in mind that for this special case the semiclassical 
approximation is exact. 

In Fig.l we first look at quantities in the time domain, i.e. as a function of n. The correlations |C„(z,z)p display 
for short times Gaussian peaks at the n-periodic points. 

However we have seen that the irreducible information about eigenstates is carried by more complicated operators 
such as the X„ in (|l^) or the Yn in (^9|). These operators carry the phase space information that corresponds to the 
pseudo-orbits in the resummation of trace formula. In fact many different powers of U contribute to them and give 
rise to complex interference structures in phase space. The pseudo-orbit contributions -contained in the coefficient 
/?„- are obtained from the fact that /?„ — {N — n) trXn ■ 

In Fig. 2 we have displayed the first few operators UXn appearing in the expansion of 5a (s)- For very short times 
we notice the simple superposition of the patterns of C„(z, z). In particular UX2 shows peaks at the period 1, period 
2 and period 3 periodic orbits. For longer times (UX^) the interferences between orbits of different periods are very 
strong and one can barely notice the period 4 orbits, while the dominant structure is of period 1. This is in striking 
contrast with the regular appearance of periodic points in the return amplitudes. Notice also the totally unexpected 
-and for us unexplained- appearance of a strong periodicity (of period three) in the phase space patterns, which is 
totally invisible in the C„ coefficients. We cannot give an explanation to this behavior except noting that a direct 
semiclassical calculation of these coefficients would be very desirable and might be more effective than one which goes 
through the calculation with periodic points. 

The smooth dependence of the operator g{s) on the quasi energy variable allow us to consider the exact response 
function as a distribution in phase space depending continuously on s, thus allowing an exact interpolation for all 
values of s. The Husimi eigen-distributions are particular values of this function which should be positive and display 
a pattern of N zero values, in accordance with the general properties of Husimi distributions of pure states on the 
torus. H 

We display in Fig. 3 the real part of 7ia(e~"^, z, z) for values of e near an eigenvalue. The distribution shows a very 
stable pattern of minima that gradually develops into TV = 14 zeros bX v = e/27r = 0.4167. At this point Z(e^") — 0, 
the distribution is positive and we verify that the imaginary part vanishes for every phase space point. This is the 
signal 1^ that a true projector has been reached and that the Husimi distribution is the modulus of an analytic 
function with N zeros. Past the eigenvalue the distribution can become negative but the pattern of minima remains 
very stable. At the values where Z{e~^^^ has an extremum, in between zeros, we show in the appendix that the 
distribution becomes flat with value 1/7V. At this point the pattern of minima changes and picks up the structure of 
zeros of the next eigenfunction. 

For cat maps the semiclassical formulae are exact and therefore the above results merely confirm that the whole 
scheme is consistent and accurate but they do not really test the semiclassical approximation. Towards this purpose 
we consider another simple linear map -the baker's- where the semiclassics is not exact. The classical T-iteration of 
the baker map from the (goiPo) to (qtiPt) real coordinates |po| , takes the complex form 

zt = Ut zq + Vt zq — wt (42) 

Zt — Vt Zq + Ut Zq — Wt (43) 
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-T 



Ut = = cosh(Tln2) (44) 

2'^ - 2-^ 

Vt = = sinh(T In 2) (45) 



Wt = —j= (7 -I- i 2 ^7^) 
v2 



(46) 



7 is the integer part of 2"^ go a-nd 7^ is the number obtained after inversion of their binary digits. Each possible 7 
defines a periodic point z-y = (7 — 17^) /\/2{2^ — 1). The return amplitudes C„(2:, z) are expressed as a sum over them 
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C„(z,z) 



E 

7=1 



exp 



(47) 



where the complex generating function takes the standard form 



g^zrz) 



1 

UT 



Vt 



(z — z^)^ + {ut — l)(z — z^)(z — z^) 



-{z-z^'f 



(48) 



and 



7 7' 



(49) 



is the action of the corresponding periodic orbit. 

In Fig. 4 we show the exact passage of the distribution through an eigenvalue {N = 16) where - as expected - the 
same general pattern as the previous example can be observed. 

In Figs. 5 and 6 we test the semiclassical approximation by comparing both the spectrum and the eigendistributions 
with the exact results. Fig. 5 shows the Z function for a small stretch of the unit circle comprising three zeroes (N = 20) 
and showing that the semiclassical zeroes are approximated well within a fraction of the mean eigenvalue separation. 
In Fig. 6 we show the three eigendistributions: the top three are the exact ones obtained at the exact zeroes, while the 
bottom three are the corresponding semiclassical ones (computed at the semiclassical zeroes). Considering that the 
distributions are obtained by combining the classsical information of 2^/^ « 1000 periodic points, it is remarkable how 
well the approximation captures the essential structure of the eigenfunctions. One important point should be noticed: 
in the semiclassical case there is no guarantee that the distribution calculated at a zero of Z will correspond exactly 
to a one dimensional projector (i.e. to a pure state) and therefore its Husimi distribution may not have exactly iV 
zeroes and may be even negative. This is clearly the case in the first semiclassical state at i/ = 0.4734. The large white 
area has negative values - and is not contoured - but corresponds very well with the central area where the zeroes of 
the exact distribution are located. The positive top of the distribution is roughly well located but it looks like the 
whole distribution had been shifted towards negative values, while mantaining a very accurate correspondence of the 
maxima and minima. The second state at v = 0.5445 reproduces very well the top part of the distribution, while now 
the minima have become positive, but again maintaining a very accurate correspondence with the exact pattern of 
zeroes. The third state at i/ = 0.6088 reproduces in very fine detail the patterns of the exact state. Overall we have 
observed that this state of affairs is quite typical: the position of the maxima and minima of the semiclassical states 
are very well reproduced but the relative heights of peaks and valleys are not always accurate. 



V. CONCLUSIONS 



We have presented a method that can, in principle, provide a reliable way to compute single Husimi eigendistri- 
butions and which has the same advantages -and drawbacks- as the resummation methods based on the spectral 
determinant for the calculation of single eigenvalues. The final results ( p^ ) and (|T]) for the Green operators when 
approximated semiclassically give approximations to the projectors on single eigenstates as sums of Gaussian centered 
on periodic points with periods up to half the Heisenberg time. Similar expressions, with contributions not so well 
localized, can be obtained for Wigner functions or probability distributions. The usual exponential proliferation of 
periodic orbits in a chaotic system is the main obstacle to applications in the deep semiclassical regime. These are 
however not limited to simple maps as the method has been applied also to the Bunimovich stadium in pl . 

We acknowledge many important discussions with F. Simonotti and P. Lebocuf, and partial support from ANPCyT 
program PICT97-01015 and ECOS-Secyt A98E03. 

APPENDIX A: NORMALIZED GREEN OPERATOR AT Z-FUNCTION EXTREMES 

To calculate the denominator of (ja we use a well known property of the parametrized matrices 



— det 

as 



= dot 



M{s) 



tr 



d 
ds 



M{s) M-^{s) 



(Al) 



If we use this property for M(s) = / — sU, then its determinant is the characteristic polynomial, and we have 
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— P(s) = P{s) tr [-U{i - sU)-'^ 
Then if P{s) is introduced into the trace, we obtain the transpose cofactor matrix, and we demonstrate 



tr 



u a 



that it is exactly the denominator of Qa- RecalUng the definition of Z{s) function 

we calculate its derivative respect to the real parameter e that rounds the unit circle (s = e~'^) 



(A2) 



(A3) 



(A4) 



dZ dZ ds i 



de ds de /?Ar/2 

Then at Sa = e~'^" where Z, as a function of e, reach a maximum or a minimum, we have the following identity 



(A5) 



(A6) 



The denominator of the cja normalized Green operator is equal to the derivative of the spectral determinant. Then, 
in a diagonal basis of U, Qa specialized at a extreme of the Z function is 



5a(Sa) 



Introducing the spectral determinant of the denominator into the sum we obtain 



N 



ga{Sa) 



N 



n=l 



\^n){lpn\ ■ 



(A7) 



(A8) 



We can separate the real and the imaginary part of the prefactor, knowing that both Sa and s„ are on the unit 
circle. 



1 - S^/Sn ~~2^ 2 [1 - 3?(5„/5„)] 



1 i 3(sa/s„) 



(A9) 



Contrary to the imaginary part, we see that the real part is independent of the index n or a. Then the normalized 
Green operator at the extremes of the Z function can be written as 



n=l ^ 



\'4'n){lpn\ 



^ N ^ N 



N 



(AlO) 



(All) 



n=l 



and splits into an hermitian and antihermitian part. We have demonstrated here that the hermitian part is exactly 
the identity operator over N. Then, the real part of the expectation value of ^(sa) will be always 1/A'', independent 
of the quantum state. 
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FIG. 1. Modulo of the coherent state return amphtude for the first six powers of the Arnold cat map propagator. The torus 
is quantized with periodic boundary conditions, and TV = 128. Linear grey scale between (white), and the maximum (black) 
of each picture. 
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FIG. 2. Modulo of the coherent state mean value for the IJ X„ operators. Same parameters and scales as Fig. 1 
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i/=0.4197 i/=0.4227 i/=0.4257 



FIG. 3. Real part of Tiais, z, z), for s = oxp(— i27ri^). Linear grey scale between (white), and 1/y/N (black). Logarithmic 
level curves at 1/N, 1/iV^, .... Arnold cat map with periodic boundary conditions, and A^" = 14. There is an eigenvalue of 

the propagator near u = 0.4167 (3rd. picture). The exact and semiclassical computation gives the same values. 




i/=0.4606 i/=0.4656 v=0A10Q 



FIG. 4. The same as Fig. 3, but for the baker map, with anti-periodic boundary conditions, and N = 16. There is an 
eigenvalue of the propagator near v = 0.4606 (4th. picture). 
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FIG. 5. Exact Z function (bold line), and its semiclassical approximation (dashed line), for the baker map, with anti-periodic 
boundary conditions, and N = 20. The zeroes correspond to the eigenvalues 
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Semiclassical 
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FIG. 6. Exact and semiclassical eigenstates in correspondence with the zeroes in Fig. 5. The gray scale and the contours 
are as in Fig.3. 
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